Direct Measurement of Effective Magnetic Diffusivity in Turbulent Flow of Liquid 

Sodium 
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The first direct measurements of effective magnetic diffusivity in turbulent flow of electro- 
conductive fluids (the so-called /3-effect) under magnetic Reynolds number Rm >> f are reported. 
The measurements are performed in a nonstationary turbulent flow of liquid sodium, generated in 
a closed toroidal channel. The peak level of the Reynolds number reached Re « 3 • 10 6 , which 
corresponds to the magnetic Reynolds number Rm « 30. The magnetic diffusivity of the liquid 
metal was determined by measuring the phase shift between the induced and the applied magnetic 
fields. The maximal deviation of magnetic diffusivity from its basic (laminar) value reaches about 
50% . 

PACS numbers: 47.65.-d, 47.27.Jv, 91.25.Cw 



Small-scale turbulence plays a crucial role in cosmic 
magnetism, providing the small-scale (turbulent) MHD- 
dynamo and contributing a lot to the dynamics of large- 
scale magnetic fields. The mean field (large-scale) dy- 
namo equations are derived by applying the Reynolds ap- 
proach to the magnetohydrodynamics (MHD) equations, 
and in the framework of the simplest case of homoge- 
neous and isotropic (but mirror asymmetric) turbulence 
they can be reduced to [l[ 
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where U and B describe the mean (large-scale) velocity 
and magnetic fields, r\ = (p/i) -1 + j3 is the magnetic dif- 
fusivity (p - electrical resistivity, fi - magnetic permeabil- 
ity), and a and /3 are the turbulent transport coefficients, 
describing the action of small-scale turbulent pulsations 
on the mean field dynamics (see e.g. @, Q)- Coefficient 
a describes the generation effects, and j3 describes the 
contribution of turbulence to diffusion of the large-scale 
magnetic field. The knowledge of the magnetic turbu- 
lent transport coefficients a and /3 is basic for astro- and 
geophysical applications in dynamo theory 

Over the last decade, major efforts were directed to- 
ward the study of MHD-dynamo in laboratory experi- 
ments (see for review Q). The first-generation dynamo 
experiments are designed on the basis of strictly-specified 
large-scale flow. The Riga dynamo is driven by the cylin- 
drical screw flow 0, the Cadarache dynamo is based on 
a von Karman flow between two counterrotating disks 
Q and even the Karlsruhe dynamo, defined as a "two- 
scale" dynamo, is driven by a set of strictly prescribed 
helical jets inside 52 tubes 0. In this sense, all labo- 
ratory dynamos can be classified as quasi-laminar. In 
spite of that, the Reynolds numbers reached about 10 7 
and the flows were fully turbulent in all experiments. 
Thus, the role of turbulence is reduced in these exper- 
iments to enhancement of the diffusion of the magnetic 
field which, with a constant magnetic permeability, can 



be considered as an increase in effective resistance of liq- 
uid metal. The growth of resistivity can be crucial for 
dynamo experiments because of the corresponding reduc- 
tion in magnetic Reynolds number. However, no direct 
measurement of effective resistivity in dynamo facilities 
has been performed up to now. An indirect indication 
of the beta-effect has been obtained in Madison sodium 
facilities by comparison of measured magnetic field and 
magnetic field simulated on the base of measured mean 
velocity field Q . An interesting scheme of eddy diffusiv- 
ity estimation from hydromagnetic Taylor-Couette flow 
experiment, recently suggested in [Toj | . 

The direct measurements of (3 are impeded by the fact 
that the effect appears only under very large Reynolds 
numbers, when numerous side-effects prevent the accu- 
rate isolation of the /3-effect. The first attempt of such 
measurements was done in a flow generated by a propeller 
in a vessel containing liquid sodium [ll|, though the au- 
thenticity of the obtained data is questionable both with 
respect to the level of the observed conductivity varia- 
tions and the estimates of the measurement errors. 

A promising method of designing high Reynolds num- 
ber flows (although nonstationary) in the limited mass of 
liquids was proposed in (12}, in which the flow was gen- 
erated by the abrupt braking of a fast-rotating toroidal 
channel. Installation of diverters in the channel made it 
possible to create a toroidal screw flow of liquid gallium, 
in which, for the first time, was observed the a-effect, 
defined by a joint action of the gradient of turbulent pul- 
sations and large-scale vorticity (l3j. The study of the 
dynamics of the nonstationary flow in a torus without di- 
verters has shown that the development of the flow in 
the channel is attended by a strong short-time burst of 
turbulent pulsations with a peak in range on the order 
of 500 - 1000Hz 14] . This burst of small-scale turbulence 
provides an opportunity to detect the increase in effective 
resistivity of liquid metal using the low frequency alter- 
nating magnetic field (~ 100 Hz). The idea of such an 
experiment has been realized in the nonstationary flow of 
liquid gallium. The toroidal channel made from textolitc 



FIG. 1: Titanium channel and thermostatic cover. 

made it possible to get magnetic Reynolds number less 
than unity flit - 

In this paper we exploit the similar experimental 
scheme using a titanium toroidal channel of larger size, 
filled with liquid sodium, which allowed us to increase the 
magnetic Reynolds number by two orders of magnitude. 

The apparatus is an electro-mechanical construction 
mounted on a rigid frame, which is used as a support 
for a rotating toroidal channel (Fig. [l| . The torus radius 
is R = 0.18 m; the radius of the channel cross-section 
is r = 0.08 m. The channel was filled with sodium in 
the vacuum and was placed into the air thermostat. The 
channel temperature may be stabilized in the range (50 — 
150)°C. The temperature sensor is mounted inside the 
channel and has good thermic contact with the sodium 
in both liquid and solid states. 

The channel is fastened on the horizontal axis, which 
is also used for mounting a driving pulley, a system of 
sliding contacts and a disk braking system. The fre- 
quency of the channel's rotation is up to 45 r.p.s. and 
the flow in the channel is generated by abrupt brak- 
ing - the braking time is no more than 0.3 sec. The 
maximum velocity of the flow is reached after channel 
is stopped and achieves about 70% of the linear veloc- 
ity of the channel before braking. This means that the 
Reynolds number Re = Ur/v [y is the kinematic viscos- 
ity of the liquid sodium) reaches at maximum the value 
Re » 3- 10 6 , which corresponds to the magnetic Reynolds 
number Rm = Urp^i w 30. 

The data-gathering system is based on an NI Data Ac- 
quisition System and is a part of an electrical measuring 
system, whose schematic circuit is shown in Fig. [5J The 
'Generator/ Amplifier' block creates in the toroidal coil 
a stabilized sinusoidal current with frequency 30 < v < 
1000 Hz, which produces an alternating toroidal mag- 
netic field inside the channel. Besides the toroidal coil, 
two diametrically located magnetic-test coils are wound 
around the channel. 



FIG. 2: Schematic circuit of the measuring system. Ud P , 
Ut, E and U are the driving pulse, the tachometer signal, 
the electromotive force of magnetic-test coil and the voltage, 
corresponding to applied current. 
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FIG. 3: Phase shift versus sodium temperature in the channel 
at frequency v = 97 Hz: experiment (points) and simulations 
(solid line). 

The change in phase shift 8 between the measured 
magnetic field and the alternating current in the toroidal 
coil is a value, which can be treated as a measure of log- 
arithmic changes of diffusivity of the sodium 

A6^C^ = C^>, (2) 
P V 

where C is a dimensional coefficient, which depends on 
the geometry and resistivity of the channel wall, and on 
the frequency of the applied magnetic field. The measur- 
ing system is completed with software, based on wavelet 
analysis, which provides calculation of the time depen- 
dence of phase shift after signals recording. Wavelets are 
required because the variation of phase shift occurs at 
times comparable with the oscillation period. 

The measurement system has been tested and cali- 
brated by measuring the dependence of the sodium re- 
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sistivity on the temperature. The channel containing the 
sodium was cooled down from 105°C to 80°C. This range 
of temperature includes the sodium freezing point, which 
gives the best measure for calibration because the resis- 
tivity of the sodium decreased at that point by 31% per- 
cent, while the temperature remained constant. This ex- 
cludes the influence of resistivity variation of titanium, 
coils, etc. Fig. [3] shows the results of phase shift mea- 
surements performed at frequency v — 97 Hz, together 
with results of numerical simulations. For this frequency, 
the skin layer thickness of titanium is about 44 mm (the 
mean thickness of titanium wall is about 10 mm) and the 
skin layer thickness of sodium is about 16 mm. 

Theoretical phase shift in the skin layer of an infinite 
cylindrical solenoid, which includes a titanium cylinder 
tube with sodium, fits the experimental points well, and 
allows us to define the factor of proportionality in relation 
J2]) for each applied frequency. For the case v = 97Hz, 
shown in Fig. 02 C — 102 ± 3 mrad. For verification 
of the method, an alternative approach of evaluating the 
sodium resistivity was used, based on the equivalent elec- 
trotechnical schematic of the transformer with the short- 
circuited secondary winding, which gave close results. 

All dynamical experiments concerning the turbulent 
flow of liquid sodium were performed under the fixed 
temperature T — (102 ± 1)°C. The estimation of sodium 
heating due to energy dissipation in decaying turbulent 
flow at the highest rotational velocity / = 50 r.p.s. (con- 
sidering that its entire kinetic energy will dissipate in the 
heat) gives AT « 0, 8°C, which corresponds to variations 
of resistivity less than 0.5%. 

Results and Summary. The rotational velocity f2 
varied from 10 to 45 r.p.s. with a step of 5 r.p.s. Mea- 
surements for all n were performed using three different 
frequencies v (53, 66 and 97 Hz). The evolution of the 
phase shift, measured at frequency v = 97 Hz for differ- 
ent velocity of channel rotation J7, is shown in Fig. |4] 
Each curve is the result of averaging over 10 realiza- 
tions. The end of braking is defined as the reference time 
point (t = 0). One can see that braking generates the 
turbulent flow, the maximal intensity of which coincides 
with the end of braking. At this moment the phase shift 
also reaches its maximum. Later on, turbulent pulsations 
rapidly decay and the phase shift reduces to zero. 

The inset of Fig. @] shows that the measured phase 
shift decays exponentially, which contradicts the ideas 
about the free decay of developed turbulence, which are 
rested on the power laws. The turbulent boundary layer 
in the nonstationary toroidal flow is developed in a very 
specific way. This was found in studies of the dynamics of 
a similar flow of liquid gallium, which have shown that 
the decay of the mean energy of the turbulent flow in 
the toroidal channel follows the t~ 2 law, while the burst 
of turbulent pulsations attends the flow formation and 
deceases abruptly 14J. This is an additional argument 




FIG. 4: Phase shift variations with flow evolution for channel 
rotation rate il = 10, 15, 40, 45 r.p.s. (bottom-up) in linear 
and lin-log (inset) scales, v = 97 Hz. 




to suppose that the measured phase shift is mostly caused 



FIG. 5: Phase shift variations with flow evolution for channel 
rotation rate = 40 r.p.s. and different frequency v, shown 
near each curve. 



by small-scale turbulence, but not by the dynamics of the 
mean flow. 

We have examined the flow across a broad range of fre- 
quencies, 31 < v < 516 Hz (the skin layer thickness varies 
then from 29 to 7 mm). Fig.[5]shows the phase shift evo- 
lution for different frequencies and it confirms the general 
idea that the turbulent diffusivity should follow the in- 
tensity of turbulent pulsation, which grows from the wall 
of the channel to its center - at a low frequency the con- 
tribution of the central part of the flow is larger and the 
/3-effect is more pronounced. At its highest frequency the 
measuring system senses only the boundary layer, which 
developed first: for v = 516 Hz the phase shift achieves 
the maximum at t = — 0.15sec, while for v = 31 Hz the 
maximum appears only at t = 0.4 sec. 

In Fig. [To] we show how the observed /3-effect depends on 
the intensity of the mean flow (on the Reynolds number, 
which is defined by the channel rotation rate before brak- 
ing). First, we show (in the upper panel) the maximal 
deviation of effective magnetic diffusivity, which corre- 
sponds to the end of braking, from the basic value. Mea- 
surements are taken using three frequencies: v — 53, 66, 
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FIG. 6: Relative increase of magnetic diffusivity (percentage) 
versus channel rotation rate fi at the end of braking (top) and 
at 0.7 sec later (bottom): v — 53 Hz (solid, black), v = 66 Hz 
(dashed, blue), v = 97 Hz (dash-dot, red). The lower panel is 
shown in logarithmic scale; dots show the power law "1/2". 



and 97 Hz. Changing frequency, we vary the depth of 
penetration of magnetic field into the turbulent flow. As 
the frequency is lowered, the thicker the skin layer be- 
comes and the more pronounced is the observed /3-effect. 
The maximal value (for f2 = 45 r.p.s. and v = 53 Hz) 
exceeds 50%. At low rotation rates the effect increases 
monotonically, in a similar manner for all frequencies; 
however, with f2 > 30 r.p.s., the monotony is disrupted 
and the curves develop in disorder. Examining individual 
curves for different realizations, it is possible to see that 
with high rotational speeds, the structure of the curve 
near the maximum becomes very complex - the maxi- 
mum becomes wider with a kind of plateau, against the 
background which appears to be separate distinct max- 
imums. All these peculiarities disappear very shortly - 
in Fig. |4] one can see that at i w 0.2 — 0.3, all curves 
evolve quite similar without any deviation. We show in 
the lower panel of Fig. [S] the deviation of effective mag- 
netic diffusivity at t = 0.7 sec. Then all three curves show 
similar monotonic increase of the /3-effect. Shown in log- 
arithmic scales, they display a tendency toward a power 
law Ar] ~ fi 1 / 2 at high rotational velocity. 

So, the measurement of electric conductivity in the 
nonstationary fully developed turbulent (Re < 3 • 10 6 ) 



flow of liquid sodium in a closed channel shows that the 
effective magnetic diffusivity essentially increases with 
the Reynolds number. For the maximal rotation rate 
Q = 45 r.p.s., which corresponds to Rm « 30, the 
maximal deviation of magnetic diffusivity reaches about 
50%. Experiments with liquid gallium at low magnetic 
Reynolds number (Rm < 1) revealed a quadratic like de- 
pendence ft ~ (Rm) 2 [Hj], which corresponds to general 
conceptions of the beta-effect for low Rm. Our results 
show that the quadratic law does not hold at moderate 
Rm. Note that the turbulent viscosity in stationary pipe 
flows at high Re increases as v t ~ Re 1,/2 [l6[ and our re- 
sults show at the highest Reynolds numbers a tendency 
to the same power law. One should treat the obtained 
dependence to the case of stationary pipe flow, or to ho- 
mogeneous turbulence, with great caution. However, in 
view of the fact that the problem of measuring the exam- 
ined characteristic in real flows is very complicated, and 
that experimental data are completely absent, measure- 
ment of the effective magnetic diffusivity in the turbulent 
medium, even in one particular case, is an important step 
toward the experimental substantiation of general MHD- 
dynamo conceptions. 

This work was supported by ISTC project 3726 and 
RFBR-SNRS grant No. 07-01-92160. 
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